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C/^ ■ Abstract 

. It has been proposed that translationally-invariant black branes are clas- 

sically stable if and only if they are locally thermodynamically stable. Reall 
' has outlined a general argument to demonstrate this, and studied in detail 

^ . the case of charged black p-branes in type II supergravity We consider 

! the application of his argument in the simplest non-trivial case, an uncharged 

' asymptotically flat brane enclosed in a finite cylindrical cavity. In this sim- 

Q ■ pie context, it is possible to give a more complete argument than in the cases 

considered earlier, and it is therefore a particularly attractive example of the 
general approach. 

1 Introduction 

It was recently proposed, in a conjecture of Gubser and Mitra P, §], that a black 
brane with a non-compact translational symmetry is classically stable if, and only if. 



X 



' it is locally thermodynamically stable. In subsequent work [|1|, Reall provided a gen- 

^ eral argument for the validity of this conjecture, based on establishing a relationship 

between the classical Gregory-Laflamme instability ^ ^ and a Euclidean negative 
mode associated with thermodynamic instability. Interesting tests of this conjecture 
are provided by solutions which pass from (classical or thermodynamic) stability to 
instability as some parameter is varied. In [|l|, charged black p-brane solutions were 
considered, where in some cases the specific heat becomes positive near extremality, 
and it was explicitly demonstrated that for suitable gauge choices, the equations satis- 
fied by the Euclidean negative mode and the classical Gregory-Laflamme perturbation 
are equivalent. 

In this paper, we will consider the simpler case of an uncharged black p-brane in d+ 
p dimensions, constructed by taking a direct product of a d-dimensional Schwarzschild 
black hole and p flat directions. We consider this system in a finite cavity, so we can 
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vary the specific heat by varying the ratio of the cavity size and the black hole's 
mass. The relation between the classically unstable mode and the Euclidean negative 
mode for this solution was already considered in infinite volume in [0]. We show that 
the introduction of the cavity walls does not spoil the equivalence, and furthermore 
show that the negative mode coincides precisely with the regime of thermodynamic 
instability.^ 

To begin, we briefly review the discussion in |T|]. In the semi-classical approxima- 
tion to the Euclidean path integral 



D[g]e~'^^^, (1) 

we write the metric as g^y = g^y + h^y, where g^y is some Euclidean solution of the 
field equations and h^y is a small perturbation. The action can then be approximated 
by 

lE[g] = Ie[9] + j d^Xy/dh'^'A^yp^hP'^. (2) 

If the operator A^ypcr in the quadratic term has negative eigenvalues, there will be an 
imaginary part in the partition function, and the classical solution gfj_y is interpreted as 
a saddle-point. If we decompose the perturbation into a transverse tracefree part hj^y, 
a trace and a longitudinal part, the trace will have a negative quadratic term, but this 
is just a sign of the usual conformal factor problem, and the resulting imaginary part 
is cancelled by a corresponding contribution from integrating over ghosts, so it does 
not correspond to a physical instability 0. The quadratic term for the longitudinal 
part is positive definite. The physical negative modes in pure gravity therefore arise 
only from the quadratic term for hj^y, which involves the Euclidean Lichnerowicz 
operator 

Gpypa = —gppgvo-^[}^^ — '^Rfipi/a- (3) 

The negative modes are given by the eigenvectors of 

G./'^/.JJ = Xhl^ (4) 

with negative eigenvalues. 

In [Q, it was shown that Schwarzschild black holes with negative specific heat 
will necessarily have a negative mode. This was extended to the charged black brane 
solutions in [|l|. Thus, thermodynamic instability in the canonical ensemble implies 
the existence of a negative mode. However, no general proof of the converse exists. 
In this paper we will check this explicitly by finding the lowest eigenvalue of and 
showing it is negative if and only if the specific heat is negative. 

The key remaining step in relating thermodynamic instability and classical insta- 
bility is relating the negative mode satisfying @ in the black hole solution to the 

similar calculation of the negative mode for Euclidean Schwarzschild-anti de Sitter was carried 
out in It was also shown in 0| that for Reissner-Nordstrom black holes in infinite volume, the 
negative mode coincides with the regime of thermodynamic instability. 
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unstable mode in the Gregory- Laflamme analysis of the perturbations about a black 
brane solution. For the uncharged black branes we are interested in, this connection 
is very direct The metric for the brane is 

ds^ = g^udx^dx'^ + 5ijdz^dz\ (5) 

where g^iy is the d- dimensional Lorentzian Schwarzschild metric and are p fiat 
spatial directions. If we consider a metric perturbation 

h^i, = exp [ifiiz') Hf,^, h^i = hij = 0, (6) 

where if^i, is transverse traceless and spherically symmetric, then the equation of 
motion implies |1 

G^u'^'^Hp^ = —fi^H^^, (7) 

where G is now the Lorentzian Lichnerowicz operator for the metric Qfj,^. A time- 
independent solution of this equation with non-zero fi represents a threshold unstable 
mode, separating stable and unstable perturbations. But on time-independent per- 
turbations, (I) and (0) are equivalent. Thus the conditions for classical instability 
and existence of a Euclidean negative mode give rise to the same equation. 

The introduction of a finite boundary introduces little modification in this ar- 
gument connecting the Euclidean negative mode in the black hole solution and the 
classical Gregory-Lafiamme instability of the uncharged black brane. For the case 
of a black hole in a finite cavity, the appropriate boundary condition in both the 
Euclidean negative mode and the classical perturbation is that the induced metric on 
the boundary be unchanged. In the remainder of the paper, we will show that this 
Euclidean negative mode (and hence the classical instability) occurs precisely when 
the black brane is thermodynamically unstable. 



2 Schwarzschild in a finite cavity 

Gravitational thermodynamics in a finite cavity was considered by York [Q. He 
considered four-dimensional spacetime, with the proper area of the spherical cavity 
A = Anrl and the local temperature T at the cavity wall fixed. In the Euclidean solu- 
tion, these boundary conditions amount to fixing the induced metric on the boundary. 
If the product of cavity radius and temperature was sufficiently high, r^T > v^/Svr, 
there were two Schwarzschild solutions of mass Mi, M2 which satisfied the boundary 
conditions. As r^T runs from this minimum value to infinity, Mi/vf, runs from 1/3 to 
0, and M2/rfe runs from 1/3 to 1/2 (thus, in the limit, the black hole fills the cavity). 
The stability in the canonical ensemble was analysed by calculating the specific heat 
at constant area of the boundary, with the result 

C,.8.M^(l--)(--l) . (8) 
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Thus, the smaller black hole of mass Mi has negative specific heat and is thermody- 
namically unstable, while the larger black hole of mass M2 has positive specific heat 
and is thermodynamically stable. 

In considering perturbations, it is more convenient to use the black hole's mass M 
and Tfo as parameters, rather than T and r^,, and we will imagine varying rj,. We can 
consider any > 2M. For 2M < < 3M, we have a positive specific heat, while for 
> 3M, we have a negative specific heat. We therefore expect to see a negative mode 
only for Vh > 3M. Allen analysed the problem of finding the negative mode for black 
holes in a finite cavity prior to York's study of the thermodynamics |T^. He found that 



the negative mode exists if ^ 2.89M. A resolution of this apparent contradiction 
was presented by York , who observed that Allen had fixed the coordinate location 
of the boundary, which, in his choice of gauge for the perturbation, did not correspond 
to fixing the area of the boundary. 

In the next section, we will discuss the solution of the eigenvalue equation (H) with 
the induced metric on the boundary fixed. First, however, we will briefly discuss the 
straightforward extension of York's thermodynamic analysis to higher dimensions. 
The (i-dimensional Schwarzschild black hole is 

, , / ujM\ , , / uM\~^ ,99,^ 



where u = IGn/Kd — 2)V^-2] and Vd-2 is the volume of the unit S'^ ^. Hence, the 
relation between mass and temperature is 

1 ujM^ "^/^ 



^^^^^ = 4vr(c.M)V(.-3) ['-4^^ ' 



and the specific heat of a black hole in a cavity of radius r?, is 

We therefore expect there to be a negative mode in the Euclidean solution only for 



3 Finding the Negative Mode 

We will now discuss the calculation of the negative mode for a d-dimensional Euclidean 
Schwarzschild black hole in a finite cavity with the induced metric on the wall of the 
cavity fixed. By the foregoing discussion, this also gives a threshold unstable mode 
for the uncharged black brane solution. We will first discuss the analysis for four 
dimensions in detail and then sketch the extension to higher dimensions. 
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3.1 Four dimensions 



We wish to consider perturbations of the four-dimensional Euchdean Schwarzschild 
metric, 



^ 2M\ ,2 dr^ 2 ,.^2 2-2 

1 1 dr^ + TTTTT + r^<i-Q + sm^ 



r 



(12) 



We are looking for perturbations which are eigenvectors of (^) with negative eigen- 
values. This problem was analysed in [§] for the black hole in infinite volume, and it 
was found that in an expansion in spherical harmonics, the only part which can have 
a negative eigenvalue is the spherically symmetric static mode. In considering the 
black hole in a finite box, it will therefore suffice for us to consider a perturbation 



= Diag 



(13) 



where the condition of transversality, V ^h^v = 0, implies 

r{r-2M) d 3r - 5M' 



Ho{r) 



3M dr 



3M 



Hiir). 



(14) 



The task of finding the negative mode then reduces to finding negative eigenvalues of 



2M 



d^ 

dr"^ 



2(r-4M)(2r-3M) d 



+ 



8M 



(r - 3M) dr r^{r - 3M) 



^i(r) 



\H^{r). 

(15) 



This equation has regular singular points at r = 0, r = 2M & r = 3M and an 
irregular singular point at r = oo, so an explicit solution to it can not be found. As 
in previous works, we must therefore use numerical methods to find the eigenvalues. 
In the infinite cavity the required boundary conditions for the solution are regularity 
at the horizon, r = 2M, and normalizability at infinity. Imposing these conditions 
on the solution then yields one negative eigenvalue, A ~ — 0.19M~^ 0. 

To search for negative modes in a finite cavity, we need to impose the condition 
that the induced metric on the walls of the cavity is fixed. This is complicated by the 
fact that the perturbation is generically non-zero in the spherical directions. In the 
unperturbed metric the area of a spherical cavity of radius r = r^, is A = Ann^. The 
perturbed metric is 



ds' 



1 - 



2M 



'1 



eHo{r))dT'+^-±^^dr' 



I — 2M 



1 



eiHoir) + H,{r)) 



and so after the perturbation the surface r = has area 

1 



I - -eiHoir) + H,ir)) 



dn^ 

(16) 



(17) 
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The change of area could be corrected by allowing the boundary of the cavity to 
move to r'f^ after the perturbation, where is chosen such that the area of the cavity 
wall in the perturbed metric is the same as that in the unperturbed metric at r^. We 
will instead give a more careful analysis by making a change of coordinates to work 
in a different gauge for the perturbation, in which the area of the cavity wall at fixed 
radius is unchanged by the perturbation. The appropriate coordinate transformation 
is 



P ■ P 



I - -e{Ho{r) + H,{r)) 



Rewriting the perturbed metric in the new coordinates and keeping only the dominant 
first order terms of the metric perturbation gives 

ds^=il- —] (1 + eFo(p)) dr^ —) (1 + eF,{p)) dp^ + p^dQ\ (19) 



where 



Mp) 

Flip) 



P 



{2p-3M)Ho{p) + MH,{p) 



P 



2 (p - 2M) 

(p - 3M) Hoip) + (3p - 7M) H,ip) + pip - 2M) (Ho'ip) + H,\p)) 



2 (p - 2M) 

Thus the metric perturbation in this gauge is 

h% = Dtag [Fo(p),Fi(p),0,0]. 



(20) 



With this new gauge we see that to hold the surface area fixed when the metric 
perturbation is switched on is a trivial task: we simply hold pb fixed. 

In the introduction, we saw that physical negative modes are determined by 
the eigenvalue equation (^), which involves only the transverse traceless part of 
the perturbation. An infinitesimal coordinate transformation changes the metric by 
9^lu 9^iv + {^liv) for some vector field S,^, and thus only changes the longitudinal 
part. That is, although the metric perturbation in our new gauge choice looks more 
complicated, the equation we need to solve is unchanged, as the transverse traceless 
part is not affected by the gauge transformation. The transverse traceless part of 
is just ([l^) with r — i> p, that is. 



h^^f"^ = Diag 



Hoip), H,{p), -1 (Hoip) + H,{p)) , -1 (Hoip) + H,ip)) 



(21) 



with the transversality condition 

p{p 



Hoip) 



2M) d 3p-5M 



P-3M dp P-3M 



^i(p)- 



(22) 
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We are thus seeking negative eigenvalues of the equation 

2M\ <f 2{p-AM){2p-3M)d^ 8M 



p J dp^ p'^{p-3M) dp p2(p-3M) 



H,{p) = \H,{p), 

(23) 



subject to the boundary conditions of regularity at p = 2M and the isothermal 
boundary condition at the fixed radial position p = p^. The first boundary condition is 
imposed by finding a regular series solution to the differential equation about p = 2M. 
There is only one such solution: 

oo 

H,{p) = J2cin{p/M-2r, (24) 



n=0 



where 



ai = -{XM^ + 2)ao, 



a2 = hxM^ + 2){2XM^ + 7)ao, 
o 

as = - — (AM2 + 5)(2A2M^ + 10AM2 + 15)ao, 
36 

a4 = — (585 + 520AM2 + 168A2M^ + 30A^M^ + 2A^M*)ao. 
360 

The isothermal boundary condition fixes the proper length around the in the r 
direction, which is given by A/^^J^OfejAr. It therefore imposes 

2M\ / 2M\ , ^ , 
1 = 1 1 + Fop, 

pb J V pb J 

^Foipk) = 

^HoiPb) = —^—H,{p,). (25) 
3M — 2p;, 

Hq is given in terms of Hi by (p2D , so the isothermal boundary condition reduces to 
a mixed boundary condition for Hi, 

Hi'jpb) 6pb'-20Mp, + 18M' ^ 

Hiipb) pk{pb-2M){2p,-3M) ■ ^ ^ 

This can be contrasted to the condition derived by Allen by holding r{, fixed WU 



Hi'in) 3n - 5M ^ 

Hlin) ^ n{n-2M) ' ^ ^ 

It should be noted that in deriving these conditions we have multiplied through by 
Pb — 3M in ([2^ ) [rf, — 3M in ([T^)]. This is acceptable so long as pb ^ 3M, but for 
Pb = 3M, we need to consider more carefully the true boundary condition 
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The method used in solving the differential equation is to numerically integrate 
from p = 2M to p = pf, using the form of the solution given by (0) as the initial data. 
However, the numerical method breaks down at p = 3M due to the singularity here 
in the differential equation. To overcome this problem we find a power series solution 
about p = 3M and find that the solution is in fact well behaved there. To evolve 
our solution through p = 3M it is therefore necessary to numerically integrate up to 
3M — 6 for some small 6 and use the results of this numerical integration to fit the 
power series solution at 3M to the data. This power series can then be evaluated at 
3M + 6 and the value of Hi and its first derivative can be extracted there to provide 
new initial conditions for a numerical solution beginning at 3M + 6 which can then 
be evolved up to pb- 

Unlike at p = 2M where only one of the series solutions was well behaved, there 
are two independent well behaved solutions at p = 3M. One of these is of order 1 
whilst the other is 0((p — 3M)^). Since the numerical solution breaks down in a 
fairly large region around 3M, we can not approach the singularity with too small a 
6. It is therefore necessary to go to fifth order in the series to maintain accuracy with 
a 5 of order 10~^ and so the particular solution to the differential equation which we 
require here is a linear combination of both series. We therefore use the values of 
Hi and Hi at 3M — 6 to find the coefficients bo and cq in the general form for the 
solution at p = 3M: 

oo 

Hlip) = 5Z + Cn{p/M - 3)=^) (p/M - 3)", (28) 

n=0 

where 

4 

bi = -gfeo, 

b.^[y-xM^)bo, 

/130 15,,, 2 9,2,.4\, 

65 = _ f 1^ + Z^AM^ - SX'M") k 

The eigenvalue spectrum is now found by a shooting method. An arbitrary value of 
A is input, the solution at p = p?, is found and tested to see if it obeys the correct 
boundary condition. If it doesn't, A is adjusted appropriately and the process is 
repeated until A is found to the required precision. We have repeated this process 
with varying cavity sizes so that the value of the lowest eigenvalue can be plotted 
against pb- 

In figure |], we give the results obtained for the boundary condition (pT]), demon- 
strating Allen's finding that the critical value of the radius is at Vb ~ 2.89M. In figure 
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Cl = — —Co, 

D 
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0, we give the results for our boundary condition (p6D, showing clearly the expected 
result that there is a negative mode only for pf, > 3M. 




3.2 Higher dimensions 

We will now briefly discuss the extension of the results above to uncharged p-branes 
in d+p dimensions for d > 4. We are interested in studying the negative modes of the 
Euclidean black hole geometry (^. In d > 4, the spherically symmetric transverse 
tracefree perturbation is 



h'^i, = Diag 



Hoir), H,{r), iH,{r) + H,ir)) , . . . , (ifo(r) + ifi(r)) 



(29) 

and the transversality condition V^/i(^ = gives 

^o{r)- 2r^-^-{d-l)u;M^'^'^ 2r'^-^ - {d - l)uM ^'^^^^ ^^^^ 
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The Euchdean negative mode is still given by the eigenvectors of (El), which becomes 



- 1 



ujM 



H'lir) 



[2d: 



,2d-6 



- wMr°'-3(3rf2 _ iirf + 18) + d{d - l)uj^M^ 



2d{d~^fuoM 



r2(2r<^-z -(^d- 1)ujM) 



r'^-2(2r'^-3 - {d-l)uM) 
H,(r) = XHUr). 



-H[{r) 
(31) 



As in the four-dimensional case, it is convenient to make a coordinate transformation 
before applying the boundary conditions at the shell. Here the appropriate transfor- 
mation is 



P ■ P 



1 - 



d-2 



(32) 



The condition that the induced metric on the shell is fixed then implies that the shell 
is at some fixed p = pb, and that 



HoiPb) 



d-3 



uM 



pt 



d-3 



tuM 



(Hoipb) + HM) = 0. 



(33) 



Using the transversality condition (|30D , this becomes a mixed boundary condition for 



2{d -l){d- 2)^^-*^ - 2(^2 -d- 2)ujMpt-^ + {d- ifoo^M^ 



HiiPb) Pbipt' - ujM)[2{d - 2)pf-3 -{d- \)uM\ 



0. (34) 



We analyse this eigenvalue problem using the same numerical methods as previously. 
The results are displayed in figure ^. In this graph the eigenvalue has been plotted 
against the cavity radius in a scaled and shifted radial coordinate, p, given by 



p{ujM)^ - 1 



(35) 



This choice is made so that in all dimensions we see that the black hole horizon is at 
p = and the specific heat (|Tl|) is negative only when p;, > 1. It is clear that in all 
cases this corresponds exactly to the existence of the negative mode. 



4 Conclusions 

We have demonstrated that an uncharged black brane in a spherical cavity is classi- 
cally unstable if and only if it is locally thermodynamically unstable. This provides 
a particularly simple and elegant example of the general connection between thermo- 
dynamic and classical instability conjectured in 0, 0]. The key element of the proof 
was showing that as we vary the volume of the cavity, the classical instability disap- 
pears at precisely the point where the specific heat changes sign. To do this, we used 
the observation that the threshold unstable mode for the classical instability is the 
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Figure 3: A(a;M)d-3 vs pb for d = 5 (closest to axis) to d = 10 (furthest from axis) 
dimensions 

analytic continuation of the Euclidean negative mode 0. We then provided the first 
analysis of the Euclidean negative mode for Schwarzschild in a finite cavity using the 
boundary conditions appropriate to the canonical ensemble. (A similar calculation of 
the negative mode for the Euclidean Schwarzschild anti-de Sitter solution was carried 
out in 0; unfortunately, it is not straightforward to construct a corresponding black 
brane solution in that case. The issue of stability for a warped black string solution 
in anti-de Sitter space was considered in |TT|].) 



One might be surprised that the presence of a boundary affects the classical in- 
stability at all; after all, the bulk metric is unchanged, so the initial behaviour of 
a perturbation with compact support near the horizon should be unaffected by the 
introduction of the boundary. However, in general relativity, the initial data are sub- 
ject to constraints, so we are not free to specify an arbitrary initial perturbation of 
compact support, and the boundary conditions can influence the allowed possibilities 
for initial perturbations. It is hence reasonable that the boundary can turn off the 
classical instability. 
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